ISSN 1175-5326 (print edition) 



Zootaxa 4586 (1): 035-064 
https://www.mapress.eom/j/zt/ 


Copyright © 2019 Magnolia Press 


Article ZOOTAXA 

ISSN 1175-5334 (online edition) 


https://doi.Org/10.11646/zootaxa.4586.l.2 

http://zoobank.Org/urn:lsid:zoobank.org:pub:2D01FAB0-53DF-4B89-A891-A51C548D4B72 


Redescription of Milnesium alpigenum Ehrenberg, 1853 (Tardigrada: Apochela) 
and a description of Milnesium inceptum sp. nov., a tardigrade laboratory model 


WITOLD MOREK 1 , ATSUSHI C. SUZUKI 2 , RALPH O. SCHILL 3 , DILIAN GEORGIEV 4 , 

MARIA YANKOVA 4 , NIGEL J. MARLEY 5 & LUKASZ MICHALCZYK 1 * 

; Institute ofZoology and Biomedical Research, Jagiellonian University, Gronostajowa 9, 30-387 Krakow, Poland 
department of Biology, Keio University School of Medicine, Hivoshi, Yokohama 223-8521, Japan 
department of Zoology. Institute of Biomaterials and Biomolecular Systems. Stuttgart University, Stuttgart, Germany 
department of Ecology and Environmental Conservation, University of Plovdiv. 24 Tzar Assen Street, 4000 Plovdiv, Bulgaria 
5 Marine Biology & Ecology Research Centre, University of Plymouth, Drake Circus, Plymouth, PL4 8AA, United Kingdom 
*Correspondence. E-mail: LM@tardigrada.net 


Abstract 

Intra- and interspecific variability, being at the very core of alpha taxonomy, has been a long-standing topic of debate 
among tardigrade taxonomists. Early studies tended to assume that tardigrades exhibit wide intraspecific variation. How¬ 
ever, with more careful morphological studies, especially those incorporating molecular tools that allow for an indepen¬ 
dent verification of species identifications based on phenotypic traits, we now recognise that ranges of tardigrade 
intraspecific variability are narrower, and that differences between species may be more subtle than previously assumed. 
The taxonomic history of the genus Milnesium, and more specifically that of the nominal species, M. tardigradum de¬ 
scribed by Doyere in 1840, is a good illustration of the evolution of views on intraspecific variability in tardigrades. The 
assumption of wide intraspecific variability in claw morphology led Marcus (1928) to synonymise two species with dif¬ 
ferent claw configurations, M. alpigenum and M quadriftdum, with M. tardigradum. Currently claw configuration is rec¬ 
ognised as one of the key diagnostic traits in the genus Milnesium, and the two species suppressed by Marcus have recently 
been suggested to be valid. In this study, we clarify the taxonomic status of M. alpigenum, a species that for nearly a cen¬ 
tury was considered invalid. We redescribe M. alpigenum, using a population collected from the locus typicus, by the 
means of integrative taxonomy, i.e. including light microscopy, scanning electron microscopy, ontogenetic observations, 
and genetic barcoding. Moreover, the redescription of M. alpigenum allowed us to verify the uncertain taxonomic status 
of two popular laboratory models that were originally considered to be M. tardigradum ; though one was recently reiden¬ 
tified as M. cf. alpigenum. Our analysis showed that both laboratory strains, despite being morphologically and morpho¬ 
metrically nearly identical to M. alpigenum, in fact represent a new species, M. inceptum sp. nov. The two species, being 
disnguishable only by statistical morphometry and/or DNA sequences, are the first example of pseudocryptic species in 
tardigrades. 

Key words: barcoding, cryptic species, integrative taxonomy, M. tardigradum s.s., phylogeny, pseudocryptic species 


Introduction 

Tardigrades, also known as water bears, are a phylum of microscopic invertebrates that dwell in marine, freshwater 
and terrestrial ecosystems (Nelson et al. 2015). The first formal descriptions of tardigrade species were published 
in the first half of the XIX century. Among them was Milnesium tardigradum Doyere 1840, described by a French 
zoologist Louis Michel Fran?ois Doyere. The species was established as the nominal taxon for the genus 
Milnesium Doyere, 1840, family Milnesiidae Ramazzotti. 1962, and the order Apochela Schuster et al., 1980. In 
the following decades, two further species of this genus were described: M. alpigenum Ehrenberg, 1853 from 
Monte Rosa (Italy/Switzerland) and M. quadriftdum Nederstrom, 1919 from Utsjoki, Finland. Both species were 
differentiated from M. tardigradum based on claw morphology (Morek et al. 2016a), a trait currently termed the 
claw configuration (CC). However, soon after M. quadriftdum was described, Marcus (1928), in his influential 
monograph on tardigrade biology, expressed an opinion that claw morphology was only a manifestation of 
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intraspecific variability, making it unsuitable for the differentiation of Milnesium species. Therefore, he 
synonymised both M. alpigenum and M. quadrifidum with M. tardigradum. This resulted in a widespread 
conviction that M. tardigradum exhibited considerable morphological variability and throughout the following 
decades researchers used to classify any Milnesium species as M. tardigradum , no matter what geographic origin of 
specimens or morphological variation ( e.g. see Ramazzotti & Maucci 1983; Dastych 1988). 

Thus, with numerous records throughout the globe, M. tardigradum inevitably became recognised as 
cosmopolitan and for decades remained the only species in the genus Milnesium. This started to change only at the 
end of the XX century, when several new Milnesium species descriptions, although based on traits other than claw 
configuration, were published (Ramazzotti, 1962; Binda & Pilato, 1990; Maucci, 1991; Pilato & Binda, 1991). 
Later, Tumanov (2006) described additional five species and proposed to standardise some of the key 
morphometric measurements. However, the greatest increase in Milnesium species descriptions occurred after the 
redescription of M. tardigradum (Michalczyk et al. 2012a, b), when 46% of all known species of the genus were 
described within only six years (i.e. from 2012 to 2017). Based on their observations, Michalczyk et al. (2012a, b) 
assumed that CC is stable at the species level and they emphasised the value of this trait in Milnesium spp. 
differentiation. However, recently Morek et al. (2016a), using an experimental approach, discovered that some 
Milnesium species may undergo ontogenetic CC change. In fact, the latest integrative study on intraspecific 
variability of M. tardigradum by Morek et al. (2019) showed that the nominal species also exhibited developmental 
variability in CC, which had gone undetected in the redescription by Michalczyk et al. (2012a, b). Importantly, 
however, the discovery of developmental variability in CC does not undermine the taxonomic value of CC itself, 
because the developmental pattern seems to be species-specific Morek et al. (2019). Moreover, the pattern of CC 
ontogenetic variability may constitute an additional set of traits for species delimitation in the genus Milnesium. 

In other words, modern advances in Milnesium taxonomy showed that the opinion of Marcus (1928) was 
incorrect, and M. alpigenum and M. quadrifidum are good species and are now pending integrative redescriptions 
(Morek et al. 2016a). These redescriptions are particularly important because with the simplistic original 
descriptions, it is impossible to differentiate M. alpigenum and M. quadrifidum from congeners that differ solely in 
quantitative (morphometric) and/or molecular traits. This, in turn, may prevent the identification of new species 
and lead to underestimation of species diversity and overestimation of species geographic ranges. Milnesium 
quadrifidum is the only known Milnesium species with the [4-4]-[4-4] CC. Therefore, describing new species that 
exhibit this claw configuration and morphology before M. quadrifidum is redescribed involves a risk of taxonomic 
inflation (Morek et al. 2016a). On the other hand, M. alpigenum, with the [3-3]-[3-3] CC, was the first described 
member of the largest group of Milnesium species defined by the CC. Thus, the redescription should verify whether 
any of the eighteen reported species with the [3-3]-[3-3] CC and smooth cuticle require synonymising with M. 
alpigenum (Morek et al. 2016a). Moreover, the two most studied Milnesium laboratory strains, one from Japan 
(Suzuki 2003) and the other from Germany (Schill et al. 2004), fit the description ofM alpigenum. Both laboratory 
strains were originally considered to be “M tardigradum ”; though the German strain was recently tentatively re¬ 
identified as “M cf. alpigenum ” (see Michalczyk et al. 2012a, Morek et al. 2016a, and Morek et al. 2019 for 
details). Thus, it is vital to verify whether the laboratory strains represent M. alpigenum, or a new species of the [3- 
3]-[3-3] CC group (Morek et al. 2016a). These laboratory strains have been used in studies on cryptobiosis (e.g. 
Hengherr et al., 2008a, b, Hengherr et al., 2009a), astrobiology (e.g. Jonsson et al. 2016), cell biology (e.g. Beisser, 
et al. 2012, Schokraie et al. 2012, Grohme et al. 2013), physiology (e.g. Reuner et al. 2010a, Forster et al. 2012), 
developmental biology (e.g. Suzuki 2003, Suzuki 2006), experimental taxonomy (e.g. Kosztyla et al. 2016, Morek 
et al. 2016b, Stec et al. 2016), and ethology (Shcherbakov et al. 2010), thus pinpointing their identity is of great 
importance. 

In this paper, we aim to clarify the taxonomy within the genus Milnesium and to verify the taxonomic status of 
the popular laboratory models. To achieve this, we integratively analysed Milnesium individuals, with the [3-3]-[3- 
3] CC and unsculptured cuticle, collected from the M. alpigenum type locality in northern Italy, and compared them 
with the two M. cf. alpigenum laboratory cultures (from Japan and Germany). Additionally, we analysed two 
similar populations, from Switzerland and Bulgaria, to test whether these represented either M. alpigenum or the 
species used to establish the laboratory strains. In either case, these additional populations could extend the genetic 
and phenotypic variability as well as the geographic range of the species in question. 

Materials and methods 

Nomenclature. Claw configuration (abbreviated throughout the text as “CC”) is denoted according to Michalczyk 
et al. (2012a, b), i. e. as a string of bracketed numbers that represent the number of points on the secondary branches 
on external and internal claws I—III, and on anterior and posterior claws IV: formula [e-i]-[a-p]. The developmental 
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terminology follows Morek et al. (2019), i.e. immature individuals: the first instar = hatchling, and the second 
instar = juvenile; mature individuals from the third instar onwards = adults. 

Sampling and specimen isolation. Detailed collection data are provided in Table 1. All moss samples were 
collected and processed according to standard methods ( e.g. Stec et al. 2015). Tardigrades were cultured following 
the protocol in Kosztyla et al. (2016), i.e. they were fed rotifers, Lecane inermis (Bryce, 1892), and kept on plastic 
Petri dishes with scratched bottoms, at a stable temperature (8 or 16°C), and in complete darkness. Individuals 
isolated from samples and/or cultures were split into analysis groups: (i) development tracking, (ii) imaging 
and morphometry in phase contrast light microscope (PCM), (iii) imaging in scanning electron microscopy (SEM), 
and (iv) DNA sequencing; see Table 1 for details. 

Milnesium alpigenum was originally described from the Monte Rosa massif, thus we sampled this locality in 
order to find a population that could be designated as the neotype series. Since the original description is very 
limited and the type material no longer exists, any Milnesium with unsculptured cuticle and the [3-3]-[3-3] CC 
found in the vicinity of the type locality could be designated as neotype M. alpigenum. Our colleagues from the 
Adam Mickiewicz University (Poland) kindly provided us with a moss sample with a candidate population 
collected in the Monte Rosa massif, from which we isolated five live females (see Table 1 for details). The small 
number of individuals was insufficient for all the required analyses, so we cultured the five females separately (to 
control for potential multiple species exhibiting similar morphology and dwelling in a single moss cushion). After 
the isolines have perpetuated for several generations, producing sufficient numbers of individuals, we sequenced 
animals from all isolines (see below for details). DNA sequencing confirmed that all five isolines represent a single 
species, thus all individuals were pooled and used for the planned analyses (Table 1). 

Alongside the Italian M. alpigenum neotype material, we also analysed a further four Milnesium populations, 
which conformed to the original description of M. alpigenum, to test whether they represent M. alpigenum or new 
species (see Table 1 for details). These were: 

• a German laboratory strain established by R.O. Schill at the University of Tubingen in 2003 with specimens 
collected in the Bebenhausen forest and subsequently maintained at the University of Stuttgart (sample code 
DE.001; “Tubingen strain”); 

• a Japanese laboratory strain established by A.C. Suzuki at the University of Keio in 2002 with specimens 
collected in Hiyoshi (sample code JR010; “Hiyoshi H-l strain”); 

• a wild population from Switzerland (sample CH.002); 

• a wild population from Bulgaria (sample BG.058). 

Microscopy and imaging. A total of 47 individuals of the neotype M. alpigenum population (IT.057), 98 from 
the Tubingen laboratory strain (DE.001), 15 from the Japanese laboratory strain (JR 010), 9 from the Swiss population 
(CH.002), and 69 from the Bulgarian population (BG.058) were used for the PCM analysis. All were mounted in 
Hoyer’s medium on microscope slides according to the recipe and protocol described by Morek et al. (2016b). 
Photographs and measurements were taken using Nikon Eclipse 50i PCM associated with a Nikon Digital Sight DS-L2 
digital camera (termed here as “a standard class/quality PCM”) and with Olympus BX53 PCM associated with 
Olympus DP74 digital camera (“high class/quality PCM”). For structures that could not be focused in a single 
photograph, a series of up to eight pictures were taken and merged into one deep-focus image using Corel Photo-Paint 
X8. 

Additionally, 20 individuals of the neotype population of M. alpigenum (IT.057), 15 from the German 
laboratory strain (DE.001), and 18 from the Bulgarian population (BG.058) were processed for SEM imaging, 
following the protocol described in Stec et al. (2015), and then examined under high vacuum in Versa 3D 
DualBeam SEM at the ATOMIN facility of the Jagiellonian University. 

Ontogenetic variability detection. In order to test for ontogenetic variability, developmental tracking 
according to Morek et al. (2016a) was employed. In brief, exuviae with eggs were incubated individually and 
emerging hatchlings were split into three subsets: (i) mounted on permanent microscope slides in Hoyer’s medium 
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TABLE 1. Information about Milnesium populations used in the study. Bolded sample codes indicate type and neotype localities. Analyses performed: (i): culturing and development 
tracking; (ii): imaging and morphometry in PCM; (iii): imaging in SEM; and (iv): DNA sequencing. Analysis numbers correspond to those used in the text. 
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one day after hatching, (ii) reared to the juvenile stage and mounted one day after moulting, and (iii) cultured to the 
first adult instar and mounted one day after moulting. A comparison of the morphology between the first three 
instars permitted observation of any ontogenetic variability in taxonomically important traits such as the CC, 
cuticle morphology and the presence/appearance of cuticular bars under claws I—III. 

Morphometries. The number of measured specimens follow recommendations by Stec et al. (2016). 
Specimens were measured according to Tumanov (2006) and Michalczyk et al. (2012a). The pt ratio is the ratio of 
the length of a given structure to the length of the buccal tube, expressed as a percentage (Pilato 1981), in the text 
the pt values are given in italics. 

Morphometric species delineation. Despite molecular analyses indicating clear genetic differences between 
the two species addressed in this study, the morphometric ranges (both absolute and relative values) overlapped 
(see below for details). Therefore, in order to test for statistical differences between the species, we used Principle 
Component Analysis (PCA) followed by a series of Student t-tests with a-level adjusted with the Benjamini- 
Hochberg correction. The PCA allows a reduction of the original dataset components while retaining the maximum 
possible variation of data. The analysis was performed in R 3.4.2. (R Core Team 2015) by preomp function, using 
only the pt values as ratios to reduce allometric effects. The dataset had ca. 18% of missing measurements, which 
were replaced by median values of each variable in a given population to avoid losing statistical power. 
Afterwards, the variables were scaled to unit variance and zero mean to minimise the effect of different scales of 
measured traits. The results were visualised using package ggfortify (version 0.4.1, Tang et al. 2016). 

Genotyping. First, we tested whether the five isolines of the M. alpigenum type locality (IT.057) and 
populations from Germany (DE.001), Japan (JP.010), Switzerland (CH.002), and Bulgaria (BG.058) represent a 
single or multiple species. This we achieved by sequencing two variable barcodes, the nuclear Internal Transcribed 
Spacer 2 (ITS-2), and the mitochondrial Cytochrome Oxidase C subunit I (COI), for eight individuals per 
population. We established that all Italian isolines represented a single species, and that all the remaining 
populations represented another species. With this knowledge, we further sequenced two conservative nuclear 
markers, small ribosomal subunit (18S rRNA), large ribosomal subunit (28 rRNA), for four individuals from each 
species (specifically, from populations IT.057 and DE.001). 

DNA was extracted from individual tardigrades following the protocol of Chelex® 100 resin (Bio-Rad), 
extraction method by Casquet et al. (2012) with modifications by Stec et al. (2015). Primer sequences and sources 
as well as references for PCR programs are listed in Table 2. All sequences were handled in BioEditve r. 7.2.5 (Hall 
1999). COI sequences were translated into amino acids to test for potential pseudogenes. Additionally the 
uncorrected p-distances were calculated utilising MEGA 7 (Kumar et al. 2016) for comparisons between species 
and populations. All sequences were aligned using the default settings of MAFFT version 7 (Katoh et al. 2002; 
Katoh & Toh 2008). The obtained alignments were edited and checked manually in BioEdit and then trimmed to 
585 bp (ITS-2) and 509 bp (COI). 


TABLE 2. Primers and references for specific protocols for amplification of the four DNA fragments sequenced in the 
study. 


DNA 

fragment 

Primer name 

Primer 

direction 

Primer sequence (5’-3’) 

Primer source 

PCR 

programme 

18S rRNA 

SSU01_F 

forward 

AACCTGGTTGATCCTGCCAGT 

Sands et al. 

Zeller (2010) 


SSU82_R 

reverse 

TGATCCTTCTGCAGGTTCACCTAC 

(2008) 


28S rRNA 

28SF0001 

forward 

ACCCVCYNAATTTAAGCATAT 

Mironov et al. 

Mironov et al 


28SR0990 

reverse 

CCTTGGTCCGTGTTTCAAGAC 

(2012) 

(2012) 

ITS-2 

ITS2_Eutar_Ff 

forward 

GCAT C GAT GAAGAACGCAGC 

Stec et al. (2018) 

Stec et al. (2018) 


ITS2_Eutar Rr 

reverse 

TCCTCCGCTTATTGATATGC 



COI 

COIMil.tarFf 

COIMil.tarRr 

forward 

reverse 

TATTTTATTTTTGGTATTTGATGTGC 

CCTCCCCCTGCAGGATC 

Morek et al. 
(2019) 

Morek et al. 

(2019) 

Phylogenetic analysis. In order to visualise evolutionary relationships between the four Milnesium 
populations analysed in this study, phylogenetic trees using all available Milnesium ITS-2 and COI sequences were 
constructed. Thus, in addition to new sequences, the data set comprised the ITS-2 sequences for the following 
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species: M. berladnicorum Ciobanu, Zawierucha, Moglan & Kaczmarek, 2014; (KT951662 from Morek el al. 
2016a); M. dornensis Ciobanu, Roszkowska & Kaczmarek, 2015; (MG923557 from Morek el al. 2019); M. 
variefidum Morek, Gasiorek, Stec, Blagden & Michalczyk, 2016a; (KT951666-7 from Morek el al. 2016a); “M 
tardigradum ”; (GQ403681-2 from Schill, unpublished); M. tardigradum sensu stricto Doyere, 1840; (MG923551- 
5 from Morek et al. 2019). Additional COI sequences for the following species were also included: M. 
berladnicorum (KT951659 from Morek et al. 2016a); M. dornensis (MG923566 from Morek et al. 2019); M. 
variefidum (KT951663 from Morek et al. 2016a); M. tardigradum s.s. (MG923558-65 from Morek et al. 2019); M. 
cf. alpigenum (KU51342 from Kosztyla et al. 2016); “M tardigradum ” (EU244603-4 from Schill, unpublished); 
“M tardigradum ” (FJ435810 from Guil & Giribet 2012); “M tardigradum” (JX683822—5 from Vicente et al. 
2013); M. sp. (EF632553 from Sands et al., unpublished); M. sp. (KX306950 from Fox et al., unpublished); M. 
spp. (KJ857001-2, KP013598, KP013601 and KP013613 from Velasco-Castrillon et al. 2015). For the outgroup, 
sequences of Diploechiniscus oihonnae Richters, 1903 (MG063724 from Gqsiorek et al. 2017a and MG923556 
from Morek et al. 2019) were used. Concatenation was run in SequenceMatrix (Vaidya et al. 2010) and the final 
ITS-2+COI alignment was 1094 bp long. 

Using PartitionFinder version 2.1.1 (Lanfear et al. 2016), under the Akaike Information Criterion (AIC), the 
best substitution model was chosen for posterior phylogenetic analysis. As COI is a protein-coding gene, the 
alignment was divided into three data blocks representing three separated codon positions. As best-fit partitioning 
scheme, PartitionFinder suggested to retain three predefined partitions for the COI data set and four predefined 
partitions for the concatenated data set separately. As RAxML (Stamatakis 2014) allows for only a single model of 
rate heterogeneity (from the GTR family) in partitioned analyses, each data set was analysed twice: first to test ah 
possible models implemented in the program (for Bayesian Inference, BI), and then for models from the GTR 
family (for Maximum Likelihood analysis, ML). The best fit-models for four partitions for BI were: TRN+I+G for 
the first codon position, K81UF+G for the second and the third codon position, and GTR+G for the ITS-2 partition. 
For ML, the best model was GTR+G. 

BI marginal posterior probabilities were calculated using MrBayes v3.2 (Ronquist & Fluelsenbeck 2003). 
Random starting trees were used and the analysis was run for ten million generations, sampling the Markov chain 
every 1000 generations. An average standard deviation of split frequencies of <0.01 was used as a guide to confirm 
that the two independent analyses had converged. The program Tracer vl .3 (Rambaut et al. 2014) was then used to 
ensure Markov chains had reached stationarity and to determine the correct “burn-in” for the analysis, which was 
the first 10% of generations. The consensus tree was obtained after summarising the resulting topologies and 
discarding the “burn-in”. In the BI consensus tree, clades recovered with posterior probability (PP) between 0.95 
and 1.00 were considered well supported, those with PP between 0.90 and 0.94 were considered moderately 
supported and those with lower PP were considered unsupported. The consensus tree was viewed and visualised by 
FigTree v. 1.4.3, available from http://tree.bio.ed.ac.uk/software/figtree. ML topologies were constructed using 
RAxML v8.0.19 (Stamatakis 2014). The strength of support for internal nodes of ML construction was measured 
using 1000 rapid bootstrap replicates. Bootstrap (BS) support values >70% on the final tree were regarded as 
significant statistical support. 

Data deposition. Raw morphometric data are deposited in the Tardigrada Register (Michalczyk & Kaczmarek 
2013) under http://tardigrada.net/register/0056.htm [M. alpigenum), http://tardigrada.net/register/0057.htm (M 
inceptum sp. nov.) as well as in Supplementary Materials SM. 1-5. The sequences of all haplotypes were uploaded 
to GenBank (accession numbers are listed under species redescription/description and in the Appendix 1). 


Results 

Taxonomic accounts 
Phylum Tardigrada Doyere, 1840 
Class Eutardigrada Richters, 1926 
Order Apochela Schuster et al., 1980 
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Family Milnesiidae Ramazzotti, 1962 
Genus Milnesium Doyere, 1840 

Milnesium alpigenum Ehrenberg, 1853 

Fig. 1, Table 3 

M. tardigradum Marcus (1928) 

Material examined: The neotype series consisting of the neotype and 37 neoparatypes (see Table 1 and "Type 
repositories” below for details). 

Integrative redescription. Females (morphometries in Table 3): Body slightly yellowish, rather slender as 
for a Milnesium (Fig. 1A). Eyes present in live specimens, quickly dissolving after fixation in Hoyer’s medium. 
Cuticle smooth in SEM and with minute pseudopores visible on the caudo-dorsal part only under high quality PCM 
(Fig. 1C). Weakly outlined pseudoplates on the caudo-dorsal cuticle visible in some specimens only under SEM 
(Fig. IB and D). Six peribuccal papillae present, with the ventral being the smallest. Six triangular peribuccal 
lamellae of unequal size; the two lateral being slightly smaller than the dorsal and ventral lamellae, i. e. with the 
4+2 configuration (identifiable only in SEM; Fig. IF). Two lateral papillae present. Buccal tube funnel-shaped 
(Fig. IE). Claws slender, primary branches with tiny accessory points, more visible on claws IV. All secondary 
branches with three points, i.e. with the [3-3]-[3-3] CC (Fig. 1G-F1). Spurs on secondary branches long and slender, 
especially on internal and anterior claws. Cuticular bars under claws I—III present. 

Males: No males were found in the sample or culture, confirming that the neotype population is 
parthenogenetic (at least facultatively). 

Juveniles: Morphologically identical to adults, except for the lack of pseudopores. 

Hatchlings: Morphologically identical to adults, except for the lack of pseudopores and the absence of 
cuticular bars under claws I—III in the majority of examined specimens (7/8 specimens = 88%). 

Ontogenetic variability: No developmental variability in the CC. Pseudopores visible only in adults. 
Cuticular bars under claws I—III mostly absent in hatchlings but always present in juveniles and adults. 

Eggs: Oval, yellow, smooth and laid in exuviae. In the culture, up to 12 eggs were recorded in a single clutch. 

DNA markers: All sequences were of a very good quality and every marker was represented by a single 
haplotype: 18S rRNA (1054 bp, MG996146); 28S rRNA (809 bp, MH000384); ITS-2 (530 bp, MH000382); and 
COI (560 bp, MH000380). Sequences are provided in Appendix 1. 

Neotype locality: 45°58T3"N, 07°57'07"E; 1370 m asl: Italy, Monte Rosa massif, lower chair-lift station of 
Macugnaga; moss on roof. 

Etymology: Ehrenberg (1853) did not explain the choice of the species name; however, it seems reasonable to 
assume that Christian Ehrenberg named the species after the Alps, the mountain chain in which the type locality, 
the Monte Rosa massif, is located. 

Type repositories: The neotype series consist of the neotype (slide IT.057.17) and 37 “neoparatypes” 
(IT.057.01-16; 18-38). The neotype and 15 neoparatypes (IT.057.01-12; 45-47) are preserved at the Institute of 
Zoology and Biomedical Research, Jagiellonian University, Gronostajowa 9, 30-387 Krakow, Poland), further 14 
neoparatypes (IT.057.13-16; 18-26) are deposited in Department of Animal Taxonomy and Ecology, Adam 
Mickiewicz University in Poznan, Umultowska 89, 61-614 Poznan, Poland, 10 neoparatypes (IT.057.27-38) are 
stored in the Marine Biology & Ecology Research Centre, University of Plymouth, Drake Circus, Plymouth, PL4 
8AA, United Kingdom, one paratype (IT.057.47) is deposited in Natural History Museum, Cromwell Road, 
London SW7 5BD, United Kingdom, and the remaining 6 neoparatypes (IT.057.39-44) are deposited in the 
collection of Binda & Pilato, Museum of the Department of Biological, Geological and Environmental Sciences, 
Section of Animal Biology “Marcello La Greca”, University of Catania, Italy. 

Phenotypic differential diagnosis: Milnesium alpigenum has the [3-3]-[3-3] CC and “smooth” cuticle (i.e. 
cuticle smooth in SEM and with minute pseudopores visible only under high quality PCM, but with no sculpturing, 
such as reticulation, on cuticle surface). This places it in the largest group of Milnesium species that share these 
characteristics (20 species; Morek et al. 2016a; Pilato & Lisi 2016; Young et al. 2016; Pilato et al. 2016; Pilato et 
al. 2017; Schlabach et al. 2018). Nevertheless, M. alpigenum differs specifically from: 
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TABLE 3. Measurements (in pm) and the pt values of selected morphological structures of 30 specimens (neotype and neoparatypes) of Milnesium alpigemim Ehrenberg, 1853 
from type locality, Monte Rosa, Italy, mounted in Hoyer’s medium. Individuals were chosen to represent the entire body length range, with as equal representation of all 
available life stages as possible. 
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External primary branch 29 15.3 - 29.9 40.0 - 50.2 21.2 44.8 4.3 2.9 24.4 43.0 

External base + secondary branch 25 10.4 - 23.0 27.8 - 34.2 15.0 31.5 3.6 1.9 18.6 32.7 

External spur 20 4.5 - 10.7 12.7 - 17.5 7.6 15.3 1.8 1.4 8.5 15.0 

. continued on the next page 








TABLE 3. (Continued) 

CHARACTER N RANGE MEAN SD Neotype 


NO 

0 s 


p 

<N 


OO 

On 


d 


^ © 
O ON 
00 <N 


NO 

O x 

<N 

r- 


p ON N© 

CO ’—* NO 


On 

<N 


d d 


43 

o 

s 

£ 


42 

& 

cd 

T3 

S3 

O 


_o 

"c3 

S-l 

43 


M- 'St- 

<N d 


On 

<N 


p <p 

<N d 


p 

d 


n® oo p 


Ln Ln 

N d 


NO 

ON 

irr 

r- 


rp 

d 


NO 

©N 

o 

NO 


in in 
d cn 


no 

d 




NO 

©N 

■^f 


<N 





<N 

'o 

*n 

<n 

On 

On 

p 

oo 

p 

*o 

d 

d 

d 1 

d 

<N 

d 

1 d 

<N 

d 1 

d 

d 

d 

1 d 

d 

d 




> 


"-h 


<n 


*n 

<n 


*o 




p 

(N 


42 


Oh 


o (U o co 

hS hS hS ^ 

c3 

u 


^ *9 NO ^ 


p 

*n 


p On 
d On 
<N 


OO 

K 


p 

(N 


p 

d 


42 

£ 

-O 

S3 

O 


35 


X X 

w w 


d 


t\ OO 

d On 
cy-j CN) 


m m in 


_o 

ts 

u 

43 


43 

o 

s 

42 

£ 


X X 

w w 


NO 

©N 

o 


42 

C3 

T3 

S3 

O 


_o 

e3 

S-H 

43 


<N on 


p 

no 


<N 

'O 


On <N 
d cNj 
*o <n 


Cl 

K 


On t-N 

d d 

'sr- <n 


on m no o t- 

■ • NCI O x 

C— *n . CN *0 in 

CO <N ON l> CO C4 


43 

o 

s 

42 

£ 


_o 

tc 

S-H 

43 


0 

9- 


42 

J_| S-H S-H S-H 

.2 .2 .2 .2 

’£h 'C 'C ’£ 

p <D <L» <D 


43 

o 

s 

42 

£ 

*C 

Oh 

S-H 

O 


42 

£ 

Td 

S3 

o 


N® 

C ' 

NO 


<*0 On 
cnJ d 


in CO CN 


N® 

ON 


On 

On 


N® 

ON 

CO 

C"~ 


NO(NcOnOcnJON(NnO 
. . . qN • • • ©N 


§h§-2 


-g 

00 

S3 


43 

O 

§ 

S-h 

42 

S-H 

o 


O 

Oh 


MILNESIUM ALPIGENUM REDESCRIPTION 


Zootaxa 4586 (1) © 2019 Magnolia Press ■ 43 








FIGURE 1. Milnesium alpigenum Ehrenberg, 1853. A—habitus, ventral view (neotype, PCM). B —habitus, dorsal view 
(SEM). C—dorsal cuticle, the arrow indicates area with visible pseudopores (neotype, PCM). D —dorsal cuticle with the barely 
visible outline of a pseudoplate (neoparatype, SEM). E —buccal apparatus (neotype, PCM). F —six peribuccal lamellae with 
the 4+2 configuration (neoparatype, SEM). G—claws 111 with the cuticular bar below (neoparatype, PCM). H —claws IV 
(neotype, PCM). All scale bars in pm. 
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• M. antarcticum Tumanov, 2006, only recorded from the Antarctic (Smykla el al. 2012), by the maximal length 
of the buccal tube (<68 pm in M. alpigenum vs >68 pm in M. antarcticum), a smaller buccal tube standard 
width (8.6-23.6 gm in M. alpigenum vs 25.9-31.8 gm in M. antarcticum), and by a statistically lower/?/ of the 
stylet support insertion point ( 61.1-70.3, on average 64.9 in M. alpigenum vs 70.0- 73. 7, on average 71.5 in M. 
antarcticum ; / 3g =16.708, p<0.001). 

• M. argentinum Roszkowska, Ostrowska & Kaczmarek, 2015, reported from Argentina, by the appearance of 
cuticle (faint pseudopores visible only with a high quality PCM on the caudal part of the dorsal cuticle in M 
alpigenum vs well-visible pseudopores in M. argentinum on the entire dorsum with a standard PCM), and by 
the lower pi of the primary branches IV (46.9-63.2 in M. alpigenum vs 28.4-36.4 in M. argentinum). 

• M. asiaticum Tumanov, 2006, recorded from Kirghizstan (type locality), China (Beasley & Miller 2007), 
Estonia (Zawierucha et al. 2014) and the Svalbard archipelago (Kaczmarek el al. 2012), by a statistically 
lower pt of primary branches III (39.7-51. 7, on average 45.8 in M. alpigenum vs 51.5-58.3, on average 55.3 in 
M. asiaticum\ ? 36 = 15.385, p<0.001) and by a lower pt of primary branches IV (46.9-63.2 in M. alpigenum vs 
63.9—76.0 in M. asiaticum). 

• M. barbadosense Meyer & Hinton, 2012, only reported from the type locality in Barbados, by a lower/?/ of the 
stylet support insertion point (61.1-70.3 in M. alpigenum vs 71.6-82.1 in M. barbadosense) and by a higher pt 
of the primary branches IV (46.9-63.2 in M. alpigenum vs 28.4-42.2 in M. barbadosense). 

• M. beatae Roszkowska, Ostrowska & Kaczmarek, 2015, only reported from the type locality in Argentina, by 
the appearance of cuticle (faint pseudopores visible only with a high quality PCM on the caudal part of the 
dorsal cuticle in M. alpigenum vs well-visible pseudopores in M. beatae on the entire dorsum with a standard 
PCM), by a more slender buccal tube (standard width/length ratio 21-38% inM alpigenum vs standard width/ 
length ratio 58-66% in M. beatae). 

• M. bohleberi Bartels, Nelson, Kaczmarek & Michalczyk, 2014, recorded from North Carolina and Tennessee, 
USA, by a more slender buccal tube (standard width/length ratio 21-38% in M. alpigenum vs standard width/ 
length ratio 54-64% in M. bohleberi). 

• M. brachyungue Binda & Pilato, 1990, recorded from the type locality in Chile and south Argentina 
(Roszkowska et al. 2016), by a higher pt of primary branches of all claws (36.0-63.2 in M. alpigenum vs 22.9— 

33.1 in M. brachyungue). 

• M. burgessi Schlabach, Donaldson, Hobelman, Miller & Lowman, 2018, reported from Kansas, USA, by a 
higher pt of the buccal tube standard width (21.0-38.4 in M. alpigenum vs 52.9-68.5 in M. burgessi) and by a 
lower pt of primary branches IV (46.9-63.2 in M. alpigenum vs 66.6-96.2. in M. burgessi). 

• M. dornensis Ciobanu, Roszkowska & Kaczmarek, 2015, recorded from Romania (type locality), Poland 
(Kaczmarek et al. 2018) and Tunisia (Gqsiorek et al. 2017b), by the appearance of cuticle (faint pseudopores 
visible only with a high quality PCM on the caudal part of the dorsal cuticle in M. alpigenum vs well-visible 
pseudopores in M. dornensis on the entire dorsum with a standard PCM) and by a statistically lower pt of the 
buccal tube standard width (21.0-38.4, on average 31.4 in M. alpigenum vs 37.8-51.6, on average 44.1 in M. 
dornensis', / 43 = 10.473, p<0.001). 

• M. eurystomum Maucci, 1991, recorded from Greenland (type locality), Chile and Argentina (Maucci 1996), 
and Mongolia (Kaczmarek & Michalczyk 2006), by a more slender buccal tube (standard width/length ratio 
21-38% inM alpigenum vs standard width/length ratio 62-65% inM eurystomum). 

• M longiungue Tumanov, 2006, reported from the type locality in the Himalayas (India) and China (Beasley & 
Miller 2007), by the presence of accessory points on primary branches, a lower pt of primary branches III 
(39.7—51 .7 in M. alpigenum vs 57.1-73.5 in M longiungue), and by a lower/?/ of primary branches IV (46.9— 

63.2 in M. alpigenum vs 81.8-92.4 in M. longiungue). 

• M minutum Pilato & Lisi, 2016, only reported from the type locality in Sicily, by a lower pt of the buccal tube 
standard width (21.0-38.4 inM alpigenum vs 38.6-42.4 in M. minutum). 

• M sandrae Pilato & Lisi, 2016, only reported from the type locality in Hawaii, by a higher pt of the stylet 
support insertion point (61.1—70.3 in M alpigenum vs 58.0-60.5 in M sandrae) and by a lower pt of the buccal 
tube standard width (21.0-38.4 in M. alpigenum vs 44.9-48.0 in M sandrae). 

• M shilohae Meyer, 2015, only reported from the type locality in Hawaii, by a lower pt of the stylet support 
insertion point (61.1-70.3 in M alpigenum vs 75.5-77.5 in M shilohae), a lower pt of the buccal tube standard 
width (21.0-38.4 in M alpigenum vs 47.1-55.9 in M shilohae), and by a higher pt of external spurs I—III 
(11.3-17.8 in M. alpigenum vs 1.9-7.5 in M shilohae). 
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• M. swansoni Young, Chappell, Miller & Lowman, 2016, only reported from the type locality in USA, by a 
higher number of peribuccal lamellae (six in M. alpigenum vs four in M. swansoni), a lower pt of the buccal 
tube standard width ( 21.0-38.4 in M. alpigenum vs 39.2-42.2 in M. swansoni), a lower pt of the posterior 
buccal tube width ( 23.7—39.4 in M. alpigenum vs 39.9-42.2 in M. swansoni), and by a lower pt of primary 
branches I ( 36.0—47.9 in M. alpigenum vs 48.4-53.7 in M. swansoni). It should be noted that the number of 
peribuccal lamellae in M. swansoni was identified only with the use of PCM, thus until SEM observations are 
made, the number of lamellae should be treated as a working hypothesis. 

• M. tumanovi Pilato, Sabella & Lisi, 2016, only reported from the type locality in Crimea, by a higher pt of the 
stylet support insertion point ( 61.1-70.3 in M. alpigenum specimens 383-983 pm long vs ca. 52.3 in M. 
tumanovi in a specimen 774 pm long) and by a lower pt of the buccal tube standard width ( 21.0-38.4 in M. 
alpigenum specimens 383-983 pm long vs ca. 55.1 in M. tumanovi in a specimen 774 pm long). 

• M. validum Pilato, Sabella, D’Urso & Lisi, 2017 only reported from the type locality in the Antarctic, 
according to the measurements presented in the description of M. validum all pt ranges overlap, but a 
comparison of specimens of a similar body length (414-509 pm in M. alpigenum and 424-482 pm in M. 
validum) shows thatM alpigenum has a shorter buccal tube (33.0-43.0 pm inM alpigenum vs 44.1—55.6 pm 
in M. validum), moreover the two species differ in the shape of secondary branches (slender in M. alpigenum 
vs robust in M. validum, compare Fig. 1G—H here and Fig. 6B-D in Pilato et al. 2017), and in the shape of 
spurs (of typical width in M. alpigenum vs very thin in M. validum). 

• M. inceptum sp. nov. (described below), recorded from Germany, Japan, Switzerland and Bulgaria—please see 
the section “Delineation of M. alpigenum and M. inceptum sp. nov.” below for a detailed differential diagnosis 
between these two pseudocryptic species. 

• M. zsalakoae Meyer & Hinton, 2010, recorded from Arizona and New Mexico (USA), by the presence of 
accessory points on primary branches and by a lower pt of primary branches of all claws ( 36.0-63.2 in M. 
alpigenum vs 64.4-102.9 in M. zsalakoae). 

Genotypic differential diagnosis: All sequences obtained for M. alpigenum were unique and distinct from the 

sequences deposited in GenBank. The ranges of the uncorrected p-distances between neotype M. alpigenum and 

sequences of other congeners are as follows: 

• 18S rRNA: 1.1%—3.6% (2.6% on average), with the most similar being M. inceptum sp. nov. from Europe 
(MH000383, present study) and the least similar being an undetermined species from Marion Island in the sub- 
Antarctic (EU266922, Sands et al. 2008). 

• 28S rRNA: 4.5%-8.0% (6.1% on average), with the most similar being an undetermined species from the 
USA (JX888585—7, Adams et al, unpublished) and the least similar being M. tardigradum s.s. from Poland 
(KC138809, Zawierucha, unpublished). 

• ITS-2: 20.4%-23.2% (20.2% on average), with the most similar being M. tardigradum s.s. from Germany 
(JF951049, Michalczyk et al. 2012a) and the least similar beingM tardigradum s.s. from France (MG923555, 
Morek et al. 2019). 

• COI: 14.8%—25.8% (17.4% on average), with the most similar beingM variefidum from the UK (KT951663, 
Morek et al. 2016a) and an undetermined species from the USA (KX306950, Fox et al.. unpublished), whereas 
the least similar being an undetermined species from the Antarctic (KP013598, Velasco-Castrillon et al. 2015). 


Milnesium inceptum sp. nov. 

Figure 2, Tables 4-5 

M tardigradum : Suzuki 2003, Schill at al. 2004, Suzuki 2006, Pfannkuchen et al. 2007, Schill 2007, Schill & Steinbruck 2007, 
Hengherr et al. 2008a, Hengherr et al. 2008b, Jonsson et al. 2008, Schill & Fritz 2008, Suzuki 2008, Takahashi et al. 2008, 
Forster et al. 2009, Hengherr et al. 2009a, Hengherr et al. 2009b, Neumann et al. 2009, Hengherr et al. 2010, Mali et al. 
2010, Reuner et al. 2010a, Reuner et al. 2010b, Schokraie et al. 2010, Shcherbakov et al. 2010, Grohme et al. 2011, 
Forster et al. 2011, Schokraie et al. 2011, Welnicz, et al. 2011, Beisser, et al. 2012, Forster et al. 2012, Schokraie et al. 
2012, Grohme et al. 2013, Wang et al. 2014, Jonsson et al. 2016. 

M. cf. alpigenum strain Mil.alp_DE.001: Kosztyla et al. (2016), Morek et al. (2016ab), Stec et al. (2016). 
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FIGURE 2. Milnesium inceptum sp. nov. A —habitus, ventral view (holotype, PCM). B —dorsal cuticle (holotype, German 
population, PCM). C —dorsal cuticle with faint pseudopores (specimen from Bulgaria, the white arrowhead indicates the area 
where the pseudopores are more densely arranged, PCM). D —dorsal cuticle with the barely visible outline of a pseudoplate 
(paratype, SEM). E —pseudoplate surface (paratype, SEM). F —buccal apparatus (holotype, PCM). G —six peribuccal lamellae 
with the 4+2 configuration. H —claws I with the cuticular bar below (paratype, PCM). I—claws IV (paratype, PCM). All scale 
bars in pm. 
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TABLE 4. (Continued) 

CHARACTER N RANGE MEAN SD Holotype 
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Posterior branches length ratio 23 61% - 81% - 68% 








Material examined: Type series consisting of 96 specimens (population DE.001) and additional 93 specimens (15 
from JP.010 population, 9 from population CH.002, and 69 from BG.058 population). See Table 1 and “Type 
repositories” below for details. 

Integrative description. Females: Body yellowish. Eyes present in live specimens, dissolved after fixation in 
Hoyer’s medium in 50% of specimens (remained visible in 1/30 = 3% specimens of the German type series, 15/15 
= 100% of the Japanese series, 9/9 specimens = 100% of the Swiss series, and in 14/23 specimens = 61% of the 
Bulgarian series). Cuticle with very small pseudopores (0.46 ±0.06pm, detectable only under a high quality PCM) 
in the German and the Swiss population and with slightly larger (but still small, 0.62 ±0.06 pm) pseudopores in the 
Bulgarian and the Japanese population (detectable under a standard PCM). In the German and the Swiss 
population, the cuticle on the entire body appears smooth under PCM (Fig. 2B), but under SEM a weak outline of a 
single dorsal pseudoplate is visible in some specimens in the caudal part of the body (Fig. 2D-E). In the Bulgarian 
and Japanese populations, no pseudoplates were detected either under PCM or in SEM. Six peribuccal papillae 
present, with the ventral being the smallest. Six triangular peribuccal lamellae of unequal size, with the two lateral 
being noticeable smaller than the two dorsal and the two ventral, i.e. with the 4±2 configuration (identifiable only 
in SEM; Fig. 2G). Two lateral papillae present. Buccal tube funnel-shaped (Fig. 2F). Primary branches with 
typically developed and clearly visible accessory points. All secondary branches with three points, i.e. with the 
[3-3]-[3-3] CC (Fig. 2H and I). Spurs on secondary branches of moderate length. Cuticular bars under claws I—III 
present in the majority of examined specimens (23/29 specimens = 79% in the German type population, 15/15 
specimens = 100% in the Japanese population, 7/9 specimens = 78% in the Swiss population, and in 11/16 
specimens = 69% in the Bulgarian population; Fig. 2H). 

Males: No males were found in German, Swiss, or Bulgarian populations and culturing of isolated virgin 
females confirmed that the type population is parthenogenetic. However, males were found to appear 
spontaneously in an otherwise parthenogenetic culture of the Japanese strain (Suzuki 2008). This suggests that the 
species is facultatively parthenogenetic with males appearing only occasionally. 

Juveniles: Morphologically identical to adults, except for the lack of the cuticular pseudopores. 

Hatchlings: Morphologically identical to adults, except for the lack of cuticular bars under claws I—III in the 
majority of examined hatchlings (14/15 specimens = 93%), and the absence of cuticular pseudopores. 

Ontogenetic variability: No developmental variability in the CC. Pseudopores visible only in adults. 
Cuticular bars under claws I—III mostly absent in hatchlings but usually present in juveniles and adults. 

Eggs: Oval, yellow, smooth and laid in the exuviae, up to 18 in a single clutch were found in laboratory 
culture. 

DNA markers: All sequences were of a very good quality. The 18S rRNA and 28S rRNA, sequenced only in 
the German type population, were 1070 bp (MH000383) and 817 bp (MH000385) long, respectively. In ITS-2, two 
haplotypes were found: HI was shared by the German, the Japanese and the Swiss population (528 bp long, 
MH000386), whereas H2 was found in the Bulgarian population 528 bp, MH000387). The p-distance between the 
two ITS-2 haplotypes was 0.8%. The COI marker exhibited three haplotypes: HI shared by the German and the 
Swiss population (658 bp, KU513422), H2 in the Japanese population (580 bp, MK628723), and H3 in the 
Bulgarian population (647 bp, MH000381). The p-distances between the COI haplotypes were as follows: 0.5% 
(HI vs H2 and HI vs H3), and 0.3% (H2 vs H3). Sequences with marked differences are provided in Appendix 1. 

Morphology and genetic markers: The sample size of four populations does not allow us to formulate strong 
conclusions on the relationship between genetic markers and animal morphology. Nevertheless, it should be noted 
that populations with COI HI (DE.001 and CH.002) exhibited statistically smaller pseudopores than populations 
with COI H2 (JP.010) and H3 (BG.058): 0.46 ±0.06 pm (DE.001) vs 0.62 ±0.06 pm (BG.058), t 2 =l A50, p<0.00\. 
No associations were observed between ITS-2 haplotypes and phenotypic taxonomic traits. 

Type locality: 48°33'42"N, 09°03'48"E; 377 m asl: Germany, Tubingen, Bebenhausen; forest; moss on soil. 

Etymology: The name of the new species originates from the Latin “inceptor ”, meaning “an initiator”, or “a 
pioneer”, as this species was among the very first tardigrade laboratory models. Milnesium inceptum sp. nov. has 
been used in a number of studies, including first studies on molecular mechanisms underlying cryptobiosis. 

Type repositories: The type series consist of the holotype (slide DE.001.34) and 96 paratypes representing 
hatchlings, juveniles and adult females (slides DE.001.01-33). The holotype (DE.001.34) with 14 paratypes 
(DE.001.04-07; 32-33) are preserved at the Institute of Zoology and Biomedical Research, Jagiellonian 
University, Gronostajowa 9, 30-387 Krakow, Poland; 18 paratypes (DE.001.08-13) are deposited in Department of 
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Animal Taxonomy and Ecology, Adam Mickiewicz University, Poznan, Umultowska 89, 61-614 Poznan, Poland; 
18 paratypes are deposited in the Department of Zoology, Institute of Biomaterials and Biomolecular Systems, 
Stuttgart University, Germany (DE.001.14-19), 18 paratypes (DE.001.20-25) are stored in Marine Biology & 
Ecology Research Centre, Plymouth University, Drakes Circus, Plymouth, PL4 8AA, United Kingdom, , one 
paratype (DE.001.34) is deposited in Natural History Museum, Cromwell Road, London SW7 5BD, United 
Kingdom, 18 paratypes (DE.001.26-31) are deposited in Department of Ecology and Environmental Conservation, 
Faculty of Biology, University of Plovdiv, Tzar Assen 24, BG-4000 Plovdiv, Bulgaria and the remaining 9 
(DE.001.01-03) are deposited in the collection of Binda & Pilato, Museum of the Department of Biological, 
Geological and Environmental Sciences, Section of Animal Biology “Marcello La Greca”, University of Catania, 
Italy. 

TABLE 5. Measurements (in pm) and the pt values of selected morphological structures of 75 specimens of Milnesium 
inceptum sp. nov. from the type locality in Tubingen (Germany) and the additional localities from Hiyoshi (Japan), 
Zurich (Switzerland), and Kazanlak Valley (Bulgaria) mounted in Hoyer’s medium. Individuals were chosen to represent 
the entire body length range, with as equal representation of all available life stages as possible. 


CHARACTER N RANGE MEAN SD 




jLim 

Pt 

jLim 

Pt 

jiim 

Pt 

Body length 

75 

326-998 

1136-1841 

628.4 

1493 

181 

164 

Peribuccal papillae length 

48 

4.4-13.1 

14.9-24.0 

8.2 

19.0 

2.1 

1.7 

Lateral papillae length 

69 

3.4-11.2 

11.8-21.6 

7.1 

16.5 

2.1 

2.1 

Buccal tube 








Length 

75 

25.8-56.2 

- 

41.6 


9.1 


Stylet support insertion point 

73 

17.1-36.4 

59.0-71.6 

26.9 

65.5 

5.4 

2.4 

Anterior width 

75 

8.4-23.8 

28.5-45.2 

15.5 

36.8 

4.4 

4.0 

Standard width 

73 

7.1-21.2 

23.1-41.7 

13.5 

32.2 

4.2 

4.5 

Posterior width 

75 

7.4-22.1 

25.2-42.7 

13.8 

32.6 

4.1 

4.4 

Standard width/length ratio 

73 

23%-42% 

- 

32% 


5% 


Posterior/anterior width ratio 

75 

74%—101% 

- 

87% 


7% 


Claw 1 lengths 








External primary branch 

65 

11.0-26.1 

34.2-51.4 

18.5 

43.2 

4.2 

3.7 

External base + secondary branch 

61 

8.3-19.7 

26.8-38.0 

12.8 

32.0 

3.1 

2.8 

External spur 

31 

2.1-6.8 

7.6-14.8 

4.4 

11.3 

1.3 

1.8 

External branches length ratio 

53 

67%—81% 

- 

73% 


4% 


Internal primary branch 

70 

10.0-24.8 

32.5-50.2 

17.5 

42.0 

4.2 

3.5 

Internal base + secondary branch 

57 

8.2-21.2 

24.4-38.8 

12.7 

31.5 

3.5 

3.1 

Internal spur 

43 

3.0-8.5 

10.8-16.8 

5.7 

13.8 

1.7 

1.5 

Internal branches length ratio 

54 

68%-88% 

- 

73% 


5% 


Claw 2 lengths 








External primary branch 

72 

10.8-28.4 

37.5-56.8 

19.6 

46.7 

4.8 

4.5 

External base + secondary branch 

59 

8.8-22.8 

27.4-41.8 

13.7 

33.5 

3.4 

2.9 

External spur 

36 

2.9-8.2 

9.2-16.6 

5.2 

12.9 

1.6 

2.0 

External branches length ratio 

57 

64%-83% 

- 

70% 


4% 


Internal primary branch 

70 

10.6-27.6 

38.8-55.0 

18.8 

45.5 

4.8 

4.3 

Internal base + secondary branch 

52 

8.3-21.0 

27.3-40.7 

13.3 

32.9 

3.7 

2.9 

Internal spur 

42 

3.2-20.1 

11.8-36.8 

6.8 

16.4 

2.8 

4.0 

Internal branches length ratio 

49 

5 9%—81 % 

- 

70% 


5% 
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TABLE 5. (Continued) 


CHARACTER 

N 

RANGE 

MEAN 

SD 


jum 

Pt 

jum 

Pt 

jum 

Pt 

Claw 3 lengths 








External primary branch 

69 

11.7-27.7 

38.8-54.5 

19.6 

46.9 

4.7 

4.3 

External base + secondary branch 

64 

8.5-21.9 

27.6-42.4 

13.9 

33.9 

3.3 

3.0 

External spur 

37 

3.3—8.7 

8.9-16.8 

5.1 

12.4 

1.6 

1.8 

External branches length ratio 

59 

64%-90% 

- 

71% 


6% 


Internal primary branch 

70 

10.6-27.2 

37.3-55.7 

18.7 

45.5 

4.6 

4.4 

Internal base + secondary branch 

42 

8.3-20.5 

25.9-39.5 

13.0 

33.2 

3.8 

2.9 

Internal spur 

44 

3.3-11.0 

11.7-21.4 

6.6 

16.2 

2.1 

2.4 

Internal branches length ratio 

40 

62%-86% 

- 

71% 


6% 


Claw 4 lengths 








Anterior primary branch 

67 

13.3-34.1 

44.2-65.7 

23.1 

55.5 

5.4 

4.7 

Anterior base + secondary branch 

63 

8.8-24.2 

31.7-46.2 

15.5 

37.5 

4.0 

3.2 

Anterior spur 

45 

2.6-11.1 

7.4-22.8 

6.0 

14.1 

2.3 

3.9 

Anterior branches length ratio 

59 

61%—76% 

- 

67% 


4% 


Posterior primary branch 

64 

12.1-32.7 

44.1-65.7 

22.5 

54.3 

5.5 

5.8 

Posterior base + secondary branch 

59 

8.5-22.8 

29.7-43.8 

15.2 

36.7 

4.0 

3.6 

Posterior spur 

48 

3.0-10.7 

11.6-20.6 

6.9 

16.4 

1.9 

2.4 

Posterior branches length ratio 

55 

60%—81 % 

- 

66% 


4% 



Phenotypic differential diagnosis. Milnesium inceptum sp. nov. has the [3-3]-[3-3] CC and “smooth" cuticle 

(, i.e. cuticle smooth in SEM and with minute pseudopores, but with no sculpturing, such as reticulation, on cuticle 

surface), which places it in the largest group of Milnesium species that share these characteristics (19 species). 

Nevertheless, M. inceptum sp. nov. differs from: 

• M. alpigenum Ehrenberg, 1853 only reported from the type locality in Italy—please see the section 
“Delineation of M. alpigenum and M. inceptum sp. nov.” below for a detailed differential diagnosis between 
these two pseudocryptic species. 

• M. antarcticum Tumanov, 2006, only reported from the Antarctic (Smykla el al. 2012), by the maximal length 
of the buccal tube (<57.0 pm in the new species vs >67.0 pm in M. antarcticum), by a lower buccal tube 
standard width (7.1-21.2 pm in the new species vs 25.9-31.8 pm in M. antarcticum), and by a statistically 
lower pt of the stylet support insertion point ( 59.0-71.6, on average 65.5 in the new species vs 70.0-73.7, on 
average 71.5 in M. antarcticunr, ? so =20.590, p<0.001). 

• M. argentinum Roszkowska, Ostrowska & Kaczmarek, 2015, recorded from Argentina, by the appearance of 
cuticle (faint pseudopores visible only with a high quality PCM on the caudal part of the dorsal cuticle in the 
new species vs well-visible pseudopores in M. argentinum on the entire dorsum with a standard PCM), the 
maximal length of the buccal tube (up to 57 pm in the new species vs up to 74 pm in M. argentinum), and by a 
lower/)? of the primary branches IV ( 44.1-65.7 in the new species vs 28.4-36.4 in M. argentinum). 

• M. asiaticum Tumanov, 2006, recorded from Kirghizstan (type locality), China (Beasley & Miller 2007), 
Estonia (Zawierucha et al. 2014), and the Svalbard archipelago (Kaczmarek el al. 2012), by a statistically 
lower/)? of primary branches IV ( 44.1-65. 7, on average 54.8 in the new species vs 63.9-76.0, on average 69.7 
in M. asiaticum', ? JV =26.040, p<0.001). 

• M. barbadosense Meyer & Hinton, 2012, only reported from the type locality in Barbados, by a lower pt of the 
stylet support insertion point ( 59.0-71.6 in the new species vs 71.6-82.1 in M. barbadosense), and by a higher 
pt of the primary branches IV ( 44.1-65.7 in the new species vs 28.4-42.2 in M. barbadosense). 

• M. beatae Roszkowska, Ostrowska & Kaczmarek, 2015, only reported from the type locality in Argentina, by 
the appearance of cuticle (faint pseudopores visible only with a high quality PCM on the caudal part of the 
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dorsal cuticle in the new species vs well-visible pseudopores in M. argentinum on the entire dorsum with a 
standard PCM), and by more elongated buccal tube (standard width/length ratio 23=42% in the new species vs 
standard width/length ratio 58-66% inM beatae ) 

• M. bohleberi Bartels, Nelson, Kaczmarek & Michalczyk, 2014, recorded from North Carolina and Tennessee, 
USA, by the more slender buccal tube (standard width/length ratio 23-42% in the new species vs standard 
width/length ratio 54-64% in M. bohleberi). 

• M. brachyungue Binda & Pilato, 1990, reported from the type locality in Chile and south Argentina 
(Roszkowska et al. 2016), by a higher pt of primary branches of claws I—III ( 32.5-56.8 in the new species vs 
22.9-27.1 in M. brachyungue) and by the pt of primary branches IV ( 44.1-65.7 in the new species vs 33.1 in 
M. brachyungue). 

• M. burgessi Schlabach, Donaldson, Hobelman, Miller & Lowman, 2018, recorded from Kansas, USA, by a 
higher pt of the buccal tube standard width ( 23.1-41.7 in the new species vs 52.9-68.5 in M. burgessi) and by 
the lower pt of primary branches IV ( 44.1-65.7 in the new species vs 66.6-96.2. in M. burgessi). 

• M. dornensis Ciobanu, Roszkowska & Kaczmarek, 2015, recorded from Romania (type locality), Poland 
(Kaczmarek et al. 2018) and Tunisia (Gasiorek et al. 2017b), by the appearance of cuticle (faint pseudopores 
visible only with a high quality PCM on the caudal part of the dorsal cuticle in the new species vs well-visible 
pseudopores in M. dornensis on the entire dorsum with a standard PCM), and by a statistically lower pt of 
buccal tube standard width ( 23.1-41.7, on average 32.2 in the new species vs 37.8-51.6, on average 44.1 in M. 
dornensis', t 22 = 10.686, p<0.001). 

• M. eurystomum Maucci, 1991, recorded from Greenland (type locality), Argentina and Chile (Maucci 1996), 
and Mongolia (Kaczmarek & Michalczyk 2006), by a more slender buccal tube (standard width/length ratio 
23-42% in the new species vs standard width/length ratio 62-65% in M. eurystomum). 

• M. longiungue Tumanov, 2006, reported from the Himalayas (India, type locality) and China (Beasley & 
Miller 2007), by the presence of accessory points on primary branches, a lower pt of primary branches III 
( 37.3—55.7 in the new species vs 57.1—73.5 in M. longiungue), and by the lower pt of primary branches IV 
( 44.1—65.7 in the new species vs 81.8-92.4 in M. longiungue). 

• M. minutum Pilato & Lisi, 2016, only reported from the type locality in Sicily, by a statistically lower pt of the 
buccal tube standard width (23.1-41.1, on average 32.2 in the new species vs 38.6-42.4, on average 41.1 in M. 
minutum', t 3 =1.990, p=0.002). 

• M. sandrae Pilato & Lisi, 2016, only reported from the type locality in Hawaii, by a higher pt of the stylet 
support insertion point ( 59.0-71.6, on average 65.5 in the new species vs 58.0-60.5, on average 58.9 in M. 
sandrae; t 22 = 8.506, p<0.001) and by a lower/)? of the buccal tube standard width (23.1-41.1 in the new species 
vs 44.9—48.0 in M. sandrae). 

• M. shilohae Meyer, 2015, only reported from the type locality in Hawaii, by a lower pt of the stylet support 
insertion point (59.0-71.6 in the new species vs 75.5-77.5 in M. shilohae) and by a lower pt of the buccal tube 
standard width (23.1-41.7 in the new species vs 47.1-55.9 in M. shilohae). 

• M. swansoni Young, Chappell, Miller & Lowman, 2016, only reported from the type locality in the USA, by a 
higher number of peribuccal lamellae (six in the new species vs four in M. swansoni', but note that the number 
of peribuccal lamellae in M. swansoni was determined only with a PCM) and by a lower pt of the buccal tube 
standard width (23.1-41.1, on average 32.2 in the new species vs 39.2-42.2, on average 40.3 in M. swansoni; 
? 7 ,=10.325, pO.OOl). 

• M. tumanovi Pilato, Sabella & Lisi, 2016, only reported from the type locality in Crimea, by a higher pt of the 
stylet support insertion point (59.0-71.6 in the new species specimens being 326-998 pm long vs 52.3 in M. 
tumanovi in a specimen 774 pm long) and by a lower/)? of the buccal tube standard width (23.1-41.7 in the 
new species in specimens 326-998 pm long vs 55.1 in M. tumanovi in a specimen 774 pm long). 

• M. validum Pilato, Sabella, D’Urso & Lisi, 2017, only reported from the type locality in the Antarctic; 
according to measurements presented in the description of M. validum all pt ranges overlap, but a comparison 
between specimens of similar body length (393—513 pm in the new species and 424-482 pm in M. validum) 
shows that M. inceptum sp. nov. has a shorter buccal tube (27.1-39.0 pm in the new species vs 44.1-55.6 in M. 
validum), moreover the two species differ in the shape of the secondary branches (typical in the new species vs 
robust in M. validum, compare Fig. 2H-I here and Fig. 6B-D in Pilato et al. 2017), and in the shape of spurs 
(moderate length and of normal width in the new species vs long and very thin in M. validum). 
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• M. zsalakoae Meyer & Hinton, 2010, recorded from Arizona and New Mexico (USA), by the presence of 
accessory points on primary branches, by a lower pt of primary branches I—III (32.5-56.8 in the new species vs 
64.4-88.6 in M. zsalakoae) and by a lower pt of primary branches IV ( 44.1-65. 7 in the new species vs 94.8— 
102.9 in M. zsalakoae). 

Genotypic differential diagnosis: Four sequences deposited in GenBank prior to this publication, labelled as 
“M tardigradum", in fact represent M. inception: two ITS-2 (GQ403681-2) and two COI (EU244603-4) (all by 
Schill, unpublished). The GQ403683 and EU244604 sequences originated from Germany and represent the same 
laboratory strain that was utilised herein to describe the new species. The sequences GQ403682 and EU244603 
originated from Japan and represent the Japanese strain, also used in the present study. 

The ranges of uncorrected p-distances between the new species and sequences of other congeners are as 
follows: 

• 18S rRNA: 1.1%—3.9% (2.9% on average), with the most similar being M. alpigenum, (MG996146, present 
study) and the least similar being an undetermined species from the USA (GQ925696, Chen et al. 
unpublished) as well as an undetermined species from South Georgia in the sub-Antarctic (EU266922, Sands 
et al. 2008). 

• 28S rRNA: 0.4%-8.8% (6.0% on average), with the most similar being an undetermined species from the 
USA (AY210826, Mallatt et al. unpublished) and another undetermined species also from the USA 
(JX888540-1, Adams et al. unpublished) and the least similar being an undetermined species from Spain 
(FJ435779-80, Guil & Giribet 2012). 

• ITS-2: 19.6%-22.8% (20.3% on average), with the most similar being M. tardigradum s.s. from Hungary and 
Poland (MG923553, Morek et al. 2019) and the least similar being M. tardigradum s.s. from France 
(MG923555, Morek et al. 2019). 

• COI: 17.8%—25.8% (19.7% on average), with the most similar being M. dornensis from Romania 
(MG923566, Morek et al. 2019) and an undetermined species from the USA (KX306950, Fox et al., 
unpublished) whereas the least similar being two undetermined species from the Antarctic (KP013601 and 
KP013598, Velasco-Castrillon et al. 2015). 


Delineation of M. alpigenum and M. inceptum sp. nov. 

The two species are genetically distinct but morphologically very similar, although not identical. Therefore, they 
could be classified as pseudocryptic species, i.e. species that can be differentiated morphologically but only with a 
detailed analysis; in this case—with the use of statistical testing of morphometric traits measured in a number of 
specimens, since the classical identification based on qualitative traits and non-overlapping morphometric ranges is 
not sufficient to tell the species apart. 

The PCA analysis indicated three components comprising 59.3% of the total variance in the pt ratios. The three 
factors were as follows: PCI: the pt of the primary and secondary branches of all claws (38.2%), PC2: the pt of 
external and posterior spurs (11.3%), and PC3: the pt of buccal tube widths and stylet support insertion point 
(9.8%). The relationships between the principal components are shown in Fig. 3. The two species did not differ in 
PCI and PC3 but in contrast, they differed in PC2, thus comparisons of the first three principal components did not 
result in congruent conclusions. Specifically, when PCI and PC3 were compared (Fig. 3B), ranges for the two 
species largely overlapped. In the PCI vs PC2 comparison (Fig. 3A), the overlap between species was smaller. 
Finally, when PC2 and PC3 were compared, the ranges barely overlapped (Fig. 3C). Thus, we compared only the 
traits constituting PC2 (i.e. the pt values of claw spurs) with a series of /-tests and adjusted a-levels. Student’s t- 
tests revealed significant differences in three of the eight compared traits (mean values, ±SD, and [ranges] for M. 
alpigenum vs M. inceptum sp. nov.): 

• spur on external claw I: 13.3 ±1.5 [11.3-17.8] vs 11.3 ±1.3 [7.6-14.8), / 46 =4.002, /K0.001; 

• spur on external claw II: 15.3 ±1.4 [12.7-17.5] vs 12.9 ±2.0 [9.2-16.6], / 54 =4.630,/><0.001; 

• spur on external claw III: 15.2 ±1.4 [12.4-17.8] vs 12.4 ±1.8 [8.9-16.8], / 4g =5.758,/?<0.001. 
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In other words, the analysis showed that M. alpigenum has statistically longer (relatively to the buccal tube 
length) external spurs thanM inceptum sp. nov. (compare also Fig. 1G-H and 2H I). 

In contrast to subtle morphometric differences, the two species exhibit considerable genetic distances in all 
four analysed markers. Specifically, they differ by 1.0% in 18S rRNA, 5.2% in 28S rRNA, 21.6% in ITS-2, and by 
18.1% in COI. Most importantly, however, the two species are not immediately related to each other (see Fig. 4 for 
the positions of both species on the Milnesium phylogenetic tree), which is the strongest evidence that the two 
species represent separate phylogenetic lineages. 

To conclude, differences both in phenotypic and genetic traits unequivocally show thatM inceptum sp. nov. is 
a bona species. Nevertheless, an extreme care must be taken when identifying these species using solely 
phenotypic data. 



PCI ( 35 . 8 %) pcl (35 8%) PC2 (12.9%) 


FIGURE 3. Graphs illustrating the relationships between the first three principal components revealed by the PCA for the 
single population of M. alpigenum Ehrenberg, 1 853 and the three pooled populations of M. inceptum sp. nov. 


Phylogenetic relationships 

Phylogenetic trees obtained with BI and ML methods did not exhibit the same topologies. Importantly, however, 
majority of nodes were weakly supported in the ML tree, thus we considered it uninformative and we focused on 
the BI tree, which had good statistical support (Fig. 4). The first clade encompasses M. tardigradum s.s. and a sister 
clade composed of M. berladnicorum, M. dornensis and M. variefidum. This group is in a sister relationship to the 
clade comprising M. inceptum sp. nov. and an undescribed species (“Milnesium hisatsinomorum'\ KX306950; Fox 
et al., unpublished). The abovementioned taxa are in polytomy with M. alpigenum and an unknown species 
(described as “M. cf. tardigradum”-, JX683822-5; Vicente et al. 2013). The following three undescribed species: 
EF632553 (Sands et al., unpublished), KP013613 and KJ857002 (Velasco-Castrillon et al. 2015) formed a clade 
that, together with another undescribed species KJ857001 (Velasco-Castrillon et al. 2015), were in polytomy with 
the abovementioned species. Finally, an undescribed species represented by two sequences (KPO13598 and 
KP013601; Velasco-Castrillon et al. 2015) was a sister group to all remaining Milnesium species. Thus, our 
analysis indicated that M. inceptum sp. nov. and M. alpigenum are not sister species, even though the relationships 
of M. alpigenum with other congeners were not fully resolved. 


Discussion 

Nearly a century since the synonymisation ofM alpigenum withM tardigradum by Marcus (1928), M. alpigenum 
is now redescribed and reinstated utilising the tools of integrative taxonomy. Moreover, the redescription of M. 
alpigenum allowed us to verify the taxonomic status of the commonly used German “Tubingen” and Japanese 
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“Hiyoshi H-l” Milnesium laboratory strains. Both strains were originally identified as “M tardigradum ” (Suzuki 
2003; Schill et al. 2004), but the redescription of M. tardigradum by Michalczyk el al. (2012a, b) showed that the 
original identifications were incorrect (see also Morek et al. 2019 for an updated delineation of M. tardigradum). In 
fact, Michalczyk et al. (2012a) hypothesised that the “Tubingen” strain might represent a new species. More 
recently, Morek et al. (2016a) tentatively classified the “Tubingen” strain as “M cf. alpigenum ” since the 
morphology of the strain fit the limited original description of M. alpigenum. However, the present study has 
shown unambiguously that both laboratory strains indeed represent a distinct taxon, M. inceptum sp. nov. 


• = concatenated ITS-2 + COI sequences 
o = only COI sequences 
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FIGURE 4. The positions of M. alpigenum Ehrenberg, 1853 and M. inceptum sp. nov. on the Milnesium phylogenetic tree 
based on the concatenated data set of ITS-2 and COI sequences. Branch support values are BI posterior probabilities. Species 
names in square brackets and in grey font are GenBank labels that are incorrect species identifications, uncertain 
identifications, or invalid names (correct identifications are provided in black font, before the incorrect labels). Filled circles 
represent both ITS-2 and COI sequences, empty circles indicate COI sequences only. Scale bar shows the number of 
substitutions per site. 

Interestingly, M. alpigenum and M. inceptum sp. nov. are barely distinguishable using standard morphometric 
traits, even though they are genetically distant (not even immediate kin; see Fig. 4). In fact, the two species would 
have been unrecognised using classical taxonomic methods alone and only the use of molecular markers allowed a 
post hoc identification of statistical phenotypic differences. Still, the two species differ morphometrically, which 
means that they are not truly cryptic. As the identified differences are minor, and can only be revealed by the use of 
statistical tests, M. alpigenum and M. inceptum sp. nov. should be termed “pseudocryptic”, meaning that they are 
species that exhibit minor morphometric differences, which are not detectable using classical diagnostic keys that 
rely on the presence or absence of qualitative (morphological) traits and non-overlapping ranges of quantitative 
(morphometric) traits. This example shows explicitly how important it is to base new species descriptions on a 
proper sample size and on both phenotypic and genetic traits. 

In our opinion, species described with a few measured individuals and without associated molecular data pose 
a serious threat to the development of tardigrade taxonomy. At the time of description (Ehrenberg 1853), M. 
alpigenum was easy to differentiate from the only congener, M. tardigradum, because there were no other known 
Milnesium species with the [3-3]-[3-3] CC. However, now there are 22 Milnesium species with this CC, with most 
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of them only differentiated by morphometries. At the time of Doyere and Ehrenberg, the knowledge and available 
analytical methods were, of course, limited compared to the present time. However, despite the undeniable increase 
in knowledge of tardigrade biology and progress in taxonomic methodology, the same process may occur again and 
again if new species are continually described with classical methods and low numbers of individuals. Even if a 
new species is obviously new (e.g. because it exhibits a unique qualitative trait) but is not integratively described, it 
may hinder the detection of similar species in the future. When describing a new species with a unique trait, it 
cannot be predicted as to whether this is the only such species in the world or whether it is the first representative of 
a complex of cryptic, pseudocryptic, or even only roughly similar taxa. Therefore, if the species is the first 
representative of a group of similar species, many of those species may go unnoticed for decades, until the nominal 
species is integratively redescribed. Furthermore, if individuals collected from distant locations that represent 
morphologically similar species are erroneously identified as the nominal species, the geographic range of the 
nominal species may be highly overestimated. 

The history of confusion with the taxonomic identity of the “Tubingen” Milnesium laboratory strain also 
underlines the vital importance of integrative redescriptions. As mentioned above, the taxonomic identification of 
the strain was possible only after the redescriptions of M. tardigradum and M. alpigenum. Importantly, however, it 
needs to be underlined that if M. alpigenum was redescribed classically, the Tubingen Milnesium laboratory strain 
would not have been recognised as a separate species, since statistical differences in a fraction of morphometric 
traits—even if they were detected—would not have been accepted by the taxonomic community as sufficient 
evidence to erect a new species. Recently, another popular laboratory strain, originally identified as “ Hypsibius 
dujardini (Doyere, 1840)”, was shown to represent a new species, Hypsibius exemplaris Gqsiorek et al., 2018. The 
reason for the misidentification of the Hypsibius strain was similar to M. inceptum sp. nov., i.e. the original 
description of H. dujardini was insufficient to identify other morphologically similar species such as H. exemplaris 
(Gqsiorek et al. 2018). As there are numerous species requiring integrative redescriptions (many of them being 
nominotypical taxa for genera and higher taxonomic ranks), we would like to urge taxonomists to prioritise 
redescribing existing taxa over describing new species (e.g. see Meier & Dikow 2004). Only when species with 
older, limited descriptions are integratively redescribed and new taxa are described integratively, may we hope to 
reliably estimate tardigrade diversity, species geographic distributions, and identify evolutionary mechanisms 
underlying these phenomena. 
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APPENDIX 1. Neotype and type DNA sequences 
Milnesium alpigenum Ehrenberg, 1853 
18S rRNA (1054 bp, MG996146): 

TAGATCGTATCATCCTACATGGATAACTGTGGTAATTCTAGAGCTAATACATGCAAAAAGCCGTCTGGCCTCGTGTCAGCGGCG 

CAGTTATTAGATTAAAACCAATATAGGCTTTTCGGGTCTATTAAACTTGTGATGAATCTGAATAACCGAAGCAAAGCGCATGGT 

CTCGTACCGGCGCTAGATCTTTCAAGTGTCTGATCTATCAGCTTGTCGTTAGGTTATGTTCCTAACGTGGCTTCGACGGGTAAC 

GGGGTATCAGGGTCCGATACCGGAGAGGGAGCCTGAGAAATGGCTACCACATCCAAGGAAGGCAGCAGGCGCGCAAATTACCCA 

CTCCCAGTTCGGGGAGGTAGTGACGAAAAATAACGATGCGGGAGCATAATGCTTCCCGTAATCGGAATGAGTACACTTTAAATC 

CTTTAACGAGGATCCATTGGAGGGCAAGTCTGGTGCCAGCAGCCGCGGTAATTCCAGCTCCAATAGCGTATATTAAAGTTGCTG 

CGGTTAAAAAGCTCGTAGTTGGATCTGGGTGTTTCGATGAGTGGTGCATCTATTCGATGTCTACTACTCCATCGACACCACAAG 

CCAACCATGTCCTTCTATACCCTTCACTGGGCGTAGATAATGGGCGGTTGGAACGTTTACTTTGAAAAAATTAGAGTGCTCAAA 

GCAGGCGTATGGCCTTGAATAATGGTGCATGGAATAATGGAATAGGACCTCGGTTCTATTTGTTGGTTTTCAAGAGCTCGAGGT 

AATGATAAAGAGGAACAGACGGGGGCATTCGTATTGCGACGTTAGAGGTGAAATTCTTGGATCGTCGCAAGACGAACTACTGCG 

AAAGCATTTGCCAAGAATGTTTTCATTAATCAAGAACGAAAGTTAGAGGTTCGAAGGCGATCAGATACCGCCCTAGTTCTAACC 

ATAAACGATGCCAACCAGCGATCTGTCGGTGTTTATTTAACGACTCGACAGGCAGCTTCCGGGAAACCAAAGTGCTTAGGTTCC 

GGGGGAAGTATGGTTGCAAAGCTGAAACTAAAGGAATGACGAAGCC 

28S rRNA (809 bp, MH000384): 

TACTAAGCGGGGGAAAAGAAACCAACGGGGATTCTCCTAGTAACTGCGAGTGAACGGAGAAAAGCCCAGCGCTGAATCCTGTAG 

CTGGTAACGGTTATGGGAGCTGTAGCGTGAAGAAGGTGTACAACCATTGCAGTCAATACACGTAAGTCTCCCTGAGTGGAGCTC 

CATCCCAAGGAGGGTGCAAGGCCCGTATCGTGTTTGACGCGTGATGGTATAGCATCTTCAGAGAGTCGGGTTGTTTGGGATTGC 

AACCTAAAGCCGGTGGTAAACTCCATCGAAGGCTAAATATGACCACGAGTCCGATAGCGAACAAGTACCGTGAGGGAAAATTGA 

AAAGCACTTTGAAGAGAGAGCGAAATAGTGCGTGAAACCGCTTAGAGGCAAGCAGATGGGGCCTCGAAGGTAGAGCAGCGAATT 
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CAGCTTGCATTTCTGCTAGACTACTGTCGGCGTAGAGATCGTAAGACTCTTGTCGATGTAGGGTGTATAGTGGAATGTGAGTGC 

ACTTTCGCTGTTTGTACGCCACCGCTGTTAAGTGTGCATCCGCTGTGGTCTTGCGTGAGGCCTTGAGTGGCTTGCTGCTCAAGT 

CACCTACGCTTGGCTATTATGCAGCGCGTTTGCCTATTAACTGGACAAGTCATTCCTATGCCAGCATCGCTTCGGTGGTGTGAT 

GTCGAACACTGGCGTGTTTATTGCTGCTTCGTGGCAGTTGACGTGCTTGCACGGCTTCAGCTGCTGGTGGTATACTGCGTCGGC 

TCTACAGGCATAGTGTAGATTTGGTGGCGAGTAGATGGCTGCCCATCTAACCC 

ITS-2 (530 bp, MH000382): 

CTTTATGAACGTTGTTTCTTCGAACGCAAATTGCGGCTATGGGTTGACCGTAGCCACGTCTGGTTGAGGGTCAAACGAAAAAAA 

AATGATAGCTACGTGTTTGCTATCGATTGTCTGTCATCCTATACTGGCCATCTCAGAGCCAGGCGAAGGCTGACAGATGAAGTA 

TCAACCCTTTGACGAGCGTATTCCTGGTCTGTAGCGGATCGGAAGCCTACGGGCGCGCACATATGCGTATATGTATGTACGGGC 

TGTGTGATGGCAGTAGGTTGGAGTCGCTGATAGGCTCTGTATCGCTTAGCTGGTTAATGCATGCGGCAGTTGTGCATATATAGA 

GCCGGCTACGAGCAGAGTTTTGACCAATCGAACGACAGCCCACTAAAGTACATACTCGCATACACAATGCGCTTGGCTTCGAGT 

ACAGATATCAGTACGCTGAATGGTCATAGGTGAAGCCATGTGTTTGTCGCTTGCGACATGCTCATACACATATACACTCATTAC 

GTTGACCTCAGCTCAGGCGAGATTAC 

COI (560 bp, MH000380): 

GCCTTTTTGTAGGTTCAGCATTGAGTATATTAATTCGACTTGAATTATCTCAACCTAATTCTATACTTATAAGAGAAGATATTT 

ATAATGCTTTTATTACAAGTCATGCTTTAGTAATAATTTTTTTTTTTGTTATACCAGTGTTAATTGGAGGGTTTGGTAATTGAT 

TGGTTCCTTTAATAATTAGTTCTCCGGATATAGCTTTTCCTCGTATTAATAATGTGAGATTTTGAATGTTAATTGCTTCTTTTT 

TGTTATTAATTTTAAGAATATTTTCTGGAGCTGGTGTAGGAGCAGGTTGAACTTTGTATCCTCCTCTTACTAATATTTATGGTC 

ATAGAAATTCTTCAGTTGATTTTGCGATCTTATCATTGCATATTGCTGGTGCTTCTTCTGTGTTTAGAGCAATAAATTTTTTAA 

CTACAATTTTTAATATACATTATTTTGGTCTTCGAATAGATAAATTACCTTTGTTTGTATGATCAATTTTTATTACGGCAATCT 

TGTTAGTTTTGGCTTTGCCTGTTTTGGCTGGTGCTATTACAATATTAATTTTGATC 

Milnesium inceptum sp. nov. 

18S rRNA (1070 bp, MH000383): 

ATTAATCAGGTATGGGGTMCTAGATCGTATCATCCTACATGGATAACTGTGGTAATTCTAGAGCTAATACATGCAATGAAGCTG 
TCTGGCCTTGTGTCAGCAGCGCAGTTATTAGATTAAAACCAATATAGGCTTTCGGGTCTATTAACTTGTGATGAATCTGAATAA 
CCGAAGCAAAGCGCATGGTCTCGTACCGGCGCTAGATCTTTCAAGTGTCTGATCTATCAGCTTGTCGTTAGGTTATGTTCCTAA 
CGTGGCTTCGACGGGTAACGGGGTATCAGGGTCCGATACCGGAGAGGGAGCCTGAGAAATGGCTACCACATCCAAGGAAGGCAG 
CAGGCGCGCAAATTACCCACTCCCAGTTCGGGGAGGTAGTGACGAAAAATAACGATGCGGGAGCATAATGCTTCCCGTAATCGG 
AATGAGTACACTTTAAATCCTTTAACGAGGATCCATTGGAGGGCAAGTCTGGTGCCAGCAGCCGCGGTAATTCCAGCTCCAATA 
GCGTATATYAAAGTTGCTGCGGTTAAAAAGCTCGTAGTTGGATCTGGGTGTTTCGATGAGTGGTGCATCTATTCGATGTCTACT 
ACTCCATCGACACCACAAGCCAACCATGTCCTTCTATACCCTTCACTGGGCGTAGATAATGGGCGGTTGGAACGTTTACTTTGA 
AAAAATTAGAGTGCTCAAAGCAGGCGTATGGCCTTGAATAATGGTGCATGGAATAATGGAATAGGACCTCGGTTCTATTTGTTG 
GTTTTCAAGAGCTCGAGGTAATGATAAAGAGGAACAGACGGGGGCATTCGTATTGCGACGTTAGAGGTGAAATTCTTGGATCGT 
CGCAAGACGAACTACTGCGAAAGCATTTGCCAAGAATGTTTTCATTAATCAAGAACGAAAGTTAGAGGTTCGAAGGCGATCAGA 
TACCGCCCTAGTTCTAACCATAAACGATGCCAACCAGCGATCTGTCGGTGTTTATTTAACGACTCGACAGGCAGCTTCCGGGAA 
ACCAAAGTGCTTAGGTTCCGGGGGAAGTATGGTTGCAAAGCTGAACTAAAGGAATTGACGAG 

28S rRNA (817 bp, MH000385): 

TACTAAGCGGAGGAAAAGAAACCAACGGGGATTCTCCTAGTAACTGCGAGTGAACGGAGAAAAGCCCAGCGCTGAATCCTGTAG 

CTGGTAACGGTTACGGGAGCTGTAGCGTGAAGAAGGTGTACAACCATTGCAGTCAATACACGTAAGTCTCCCTGAGTGGAGCTC 

CATCCCAAGGAGGGTGCAAGGCCCGTATCGTGTTTGACGCGTGATGGTATAGCATCTTCAGAGAGTCGGGTTGTTTGGGATTGC 

AACCTAAAGCCGGTGGTAAACTCCATCGAAGGCTAAATATGACCACGAGTCCGATAGCGAACAAGTACCGTGAGGGAAAATTGA 

AAAGCACTTTGAAGAGAGAGCGAAATAGTGCGTGAAACCGCTTAGAGGCAAGCAGATGGGGCCTCGAAGGTAGAGCAGCGAATT 

CAGCTTGCATCTCTGCTGGACTACTGTTGGCGTAGAGATCGTAAGACTCTTGTCAATGTAGGCTGTCATAGTGGAATGTGAGTG 

CACTTTCGCTGTTTGTACGCCACCGCTGATAAATGTTTCTGCATCCGCTGTGGCCTTGTGTGAGGCCTTGAGTGGCTTGCTACT 

CAAGTAACCTACACTTGGCTATATACAGCGCGTTTGCCTTTTATCTGGTCGAGTCACATCCTATGCCGGCATTGCTTTACGGCG 

ATGCAGTGTAGATACTAGCGTGTTTATTGCTACTTCGCATTGCGGTTGACGTGCTTGCACGGCTGCTGCGGCTGGTGGTATACT 

GCGTTGGATCTACTGGTATAGTGTAGATTTGGTGGCGAGTAGACGGCTGCCCATCTAACCC 

ITS-2/H1, Germany (DE.001), Japan (JP.010) and Switzerland (CH.002) (528 bp, MH000386): 

CTTTATGAACGTTGTTTCTTCGAACGCAAATTGCGGCTATGGGTTGACCGTAGCCACGTCTGGTTGAGGGTCAAACGAAAAAAT 

CACTGATAGCTACGTGTTTGCTGTCGATTGTCTGTCAGTCTATACTGGCCATTTTAGTGCCAGGTCAAGGCTGACAGATGAAGT 

TTGGACCCTATGGCTAGCGTACTTTTGGTCTGTGACGGATCGGGAGCTGGCGCGTACATGCTCATATAGAGTTATACTGGGACT 

GTGTGATCACAGTAGGTTGGAGTCGCTGATGGGCTCTACCGTTTAGCTGACAATGCATGCGGCAGTTGTGCATAACTAGTCAGC 

GACGCGTAGACGCTTGGCCAACCGAACGACAGTCCACTCACGGTATACTCTGGCATGTTCAAGCGTACGGCTACCGAGTGCAGC 

TTCCAATACGCAGTGATATAGTCATAGGTTGATAAAGCGTGTGTCGTACGCTTCATGTGCGCGACGCATACACACTCATTACGT 

TGACCTCAGCTCAGGCGAGATTAC 
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ITS-2/H2, Bulgaria (BG.058) (528 bp, MH000387): 

CTTTATGAACGTTGTTTCTTCGAACGCAAATTGCGGCTATGGGTTGACCGTAGCCACGTCTGGTTGAGGGTCAAACGAAAAAAT 

CACTGATAGCTACGTGTTTGCTATCGATTGTCTGTCAGTCTATACTGGCCATTTTAGTGCCAGGTCAAGGCTGACAGATGAAGT 

TTGGACCCTATGGCTAGCGTACTTTTGGTCTGTGACGGATCGGAAGCCGGCGCGTACATGCTCATATAGAGTTATACTGGGACT 

GTGTGATCACAGTAGGTTGGAGTCGCTGACGGGCTCTACCGTTTAGCTGACAATGCATGCGGCAGTTGTGCATAACTAGTCAGC 

GACGCGTAGACGCTTGGCCAACCGAACGACAGTCCACTCACGGTATACTCTGGCATGTTCAAGCGTACGGCTACCGAGTGCAGC 

TTCCAATACGCAGTGATATAGTCATAGGTTGATAAAGCGTGTGTCGTACGCTTCATGTGCGCGACGCATACACACTCATTACGT 

TGACCTCAGCTCAGGCGAGATTAC 


COI/H1, Germany (DE.001) and Switzerland (CH.002) (658 bp, KU513422): 

TATATTGTATTTTATTTTTGGTATTTGATGTGCTTTTGTAGGTTCAGGTTTAAGTGTGTTAATTCGTCTTGAATTATCTCAGCC 

TAACACAATATTAATAAGTGAAGATATTTATAATGCTTTTATTACAAGTCATGCTTTAGTAATGATTTTTTTTTTTGTTATACC 

TGTTTTAATTGGAGGTTTTGGAAATTGATTAGTTCCTCTTATAATTAGATCACCAGATATAGCTTTTCCTCGTATTAACAATGT 

AAGATTTTGAATATTAGTTGCTTCTTTTGGTTTGTTGCTTTTTAGAATATTTAGGGGTACAGGAGTAGGAGCTGGTTGAACACT 

TTATCCTCCGTTAACTAGGTATAATGGCCATAGCAGTCATGCTGTCGATTATGCAATTTTGTCTTTACATATTGCAGGAGCATC 

GTCAATTTTTAGTGCACTGAATTTTTTAACGACGATTATTAATATACACTATTTTGGAGTACGAATAGATAAATTACCGTTGTT 

TGTGTGATCGATTTTTATTACTGCTCTATTGTTAGTTTTGGCTTTACCAGTACTTGCTGGAGCAATTACAATATTAATTTCTGA 

TCGTAATTTTACTACTACATTTTTTGATCCGGCAGGGGGAGGAGATCCTGTTTTATTTCAACATTTATTT 

COI/H2, Japan (JP.010) (580 bp, MK628723): 

TTTTGTAGGTTCAGGTTTAAGTGTGTTAATTCGTCTTGAATTATCTCAGCCTAACACAATATTAATAAGTGAAGATATTTATAA 
TGCTTTTATTACAAGTCATGCTTTAGTAATGATTTTTTTTTTTGTTATACCTGTTTTAATTGGAGGTTTTGGAAATTGATTAGT 
TCCTCTTATAATTAGATCACCAGATATAGCTTTTCCTCGTATTAATAATGTAAGATTTTGAATATTAGTTGCTTCTTTTGGTTT 
GTTGCTTTTTAGAATATTTAGGGGTACAGGAGTAGGAGCTGGTTGAACACTTTATCCTCCGTTAACTAGGTATAATGGTCATAG 
CAGTCATGCTGTCGATTATGCAATTTTGTCTTTACATATTGCAGGAGCATCGTCAATTTTTAGTGCACTGAATTTTTTAACGAC 
GATTATTAATATACACTATTTTGGAGTACGAATAGATAAATTACCGTTGTTTGTGTGATCGATTTTTATTACTGCTCTATTGTT 
AGTTTTGGCTTTACCAGTACTTGCTGGAGCAATTACAATATTAATTTCTGATCGTAATTTCACTACTACATTTTTT 

COI/H3, Bulgaria (BG058) (647 bp, MH000381): 

TTTTGTAGGTTCAGGTTTAAGTGTGTTAATTCGTCTTGAATTATCTCAGCCTAACACAATATTAATAAGTGAAGATATTTATAA 

TGCTTTTATTACAAGTCATGCTTTAGTAATGATTTTTTTTTTTGTTATACCTGTTTTAATTGGAGGTTTTGGAAATTGACTAGT 

TCCTCTTATAATTAGATCACCAGATATAGCTTTTCCTCGTATTAATAATGTAAGATTTTGAATATTAGTTGCTTCTTTTGGTTT 

GTTGCTTTTTAGAATATTTAGGGGTACAGGAGTAGGAGCTGGTTGAACACTTTATCCTCCGTTAACTAGGTATAATGGCCATAG 

CAGTCATGCTGTCGATTATGCAATTTTGTCTTTACATATTGCAGGAGCATCGTCAATTTTTAGTGCACTGAATTTTTTAACGAC 

GATTATTAATATACACTATTTTGGAGTACGAATAGATAAATTACCGTTGTTTGTGTGATCGATTTTTATTACTGCTCTATTGTT 

AGTTTTGGCTTTACCAGTACTTGCTGGAGCAATTACAATATTAATTTCTGATCGTAATTTCACTACTACATTTTTTGATCCGGC 

AGGGGGAGGAGATCCTGTTTTATTTCAACATTTATTTTGGTTTTTTGGGCATCCAGAAG 


Supplementary Materials 

Supplementary Materials SM.01— Morphometric measurements of the neotype (Italy, IT.057) population of Milnesium 
alpigenum Ehrenberg, 1853. 

Supplementary Materials SM.02— Morphometric measurements of type population (Germany, DE.001) of Milnesium inceptum 

sp. nov. 

Supplementary Materials SM.03 —Morphometric measurements of the Japan (JP.010) population of Milnesium inceptum sp. 

nov. 

Supplementary Materials SM.04— Morphometric measurements of the Swiss (CH.002) population of Milnesium inceptum sp. 

nov. 

Supplementary Materials SM.05 —Morphometric measurements of the Bulgarian (BG.058) population of Milnesium inceptum 

sp. nov. 
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